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The electron-nuclei (hyperfine) interaction is central to spin qubits in solid state systems. It 
can be a severe decoherence source but also allows dynamic access to the nuclear spin states. We 
study a double quantum dot exposed to an on-chip single-domain nanomagnet and show that its 
inhomogeneous magnetic field crucially modifies the complex nuclear spin dynamics such that the 
Overhauser field tends to compensate external magnetic fields. This turns out to be beneficial for 
polarizing the nuclear spin ensemble. We reach a nuclear spin polarization of c± 50%, unrivaled in 
lateral dots, and explain our manipulation technique using a comprehensive rate equation model. 
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The hyperfine interaction (HFI) between few electrons 
and a bath of nuclear spins induces a complex quantum 
many-body dynamics which has been studied in a va- 
riety of systems including phosphorus donors in silicon 
PQ, nitrogen vacancy centers in diamond [2 , quantum 
Hall systems [3], and semiconductor-based quantum dots, 
both optically [4] and in transport [5j [6] . In GaAs dou- 
ble quantum dots (DQDs), each electron interacts with 
^10 6 nuclear spins, which fluctuate thermally even at 
cryogenic temperatures. Their HFI causes electron spin 
decoherence [6 -9 , but it also offers a means to control the 
nuclear spins dynamically [T0HT3] . As has been proposed 
[T4HT7] and demonstrated [T8H20] , nuclear spin manipu- 
lation facilitates nuclear state preparation, which can en- 
hance spin coherence times [21, 22 . Nuclear spin manip- 
ulation has also motivated theoretical proposals [23j [24] 
and experimental realizations [25 , 26] of nuclear spins as 
quantum memory. 

Here, we present a novel nuclear spin manipulation 
technique. We couple a DQD with a single-domain nano- 
magnet whose inhomogeneous magnetic field does not de- 
pend on external fields. This allows measurements in a 
new regime, which will increase our fundamental under- 
standing and control of the coupled electron-nuclei sys- 
tem common to a variety of platforms. As an example, 
we demonstrate unusually strong dynamic nuclear spin 
polarization (DNSP). In laterally defined DQDs, adia- 
batic pumping experiments have produced polarizations 
of 1-5 % [TTUT3] , and using electron dipole spin reso- 
nance, 16% has been realized [10]. Using a simpler 
experimental technique, we report ^ 50 % polarization, 
achieved by exploiting the benefits of our nanomagnet. 

We measure the current, /, which results from a dc 
voltage, here V = 1 mV, applied across the DQD exposed 
to the inhomogeneous magnetic field of a nanomagnet 
(see Figures []Ji,c). As detailed in Figure [jj), electrons 
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FIG. 1. (a) Scanning electron micrograph showing the Ti/Au 
gates (yellow) used to define the DQD. Electrons tunnel se- 
quentially from source (S) to drain (D) through the two QDs 
located near •. The (blue) single-domain Co nanomagnet 
(200 (width) x 50 (height) x 2000 nm 3 , 2/5 visible) gener- 
ates an inhomogeneous magnetic field, B nm . (b) The tunnel- 
ing sequence (01) (1 1) (0 2) (01). Triplet- singlet 
transitions {T_,To,T+} — > Sn require lifting of the PSB. (c) 
The relevant layers of the heterostructure, gates, and nano- 
magnet. The magnetic field lines (blue) are calculated accord- 
ing to ref. m| The black arrows indicate the local magnetic 
fields B^ m = (0,32.3, 22.4) mT and B^ m = (0, 5.2, 10.6) mT. 
The 3 mT radius of • represents the typical Overhauser field 
fluctuations. 



tunnel sequentially through the DQD via the occupation 
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cycle (01) -> (11) H> (0 2) (01), where (ran) indi- 
cates the number, m (n), of electrons in the left (right) 
dot. The transition (01) —> (11) loads one of four (11) 
states, which consist of the singlet state, Sn, and the 
three triplet states, Tn = {T_, To, T + }. Only T± have 
a nonzero spin projection along the quantization axis, 
which we choose parallel to the external magnetic field, 
B ext (along the z-axis in Figure [IJi) . The only energeti- 
cally accessible (0 2) state is the singlet state, So2- The 
corresponding eigenenergies are plotted in Figure [2^i as 
a function of the energy detuning, A, between the dia- 
batic singlet states Sn and So2- The singlet eigenstates, 
S s and 5 a , are the symmetric and antisymmetric com- 
binations of Su and So 2 mixed by the interdot tunnel 
coupling, t c . Tq is at zero energy (neglecting exchange 
coupling), while T± are split by their Zeeman energies, 
E± = ±gfji B (|B L + B R |) /2, where B L (B R ) is the local 
magnetic field in the left (right) dot, g the Lande g-factor, 
and fiB the Bohr magneton. For A > 0, the triplets are 
well separated from S s , but not from S a . 

In a homogeneous magnetic field (B L = B R ) , tran- 
sitions between triplets and singlets are forbidden by 
Pauli spin blockade [28] (PSB) (dashed arrows in Fig- 
ure [TJd). Eventually the occupation cycle stalls in one 
of the triplets resulting in / = (neglecting cotunnel- 
ing). DNSP requires 7^0, which can be induced by 
a local field difference, B L - B R , that lifts the PSB by 
coupling triplets to singlets. One way to produce in- 
homogeneous fields is to include on-chip micromagnets, 
which have been used for all-electric control of a single 
electron spin [29H5T] . However, at external fields be- 
low a few hundred mT, micromagnets form magnetic do- 
mains, which greatly reduce their fields. Here, we use a 
nanomagnet (see Figure [IJi), which forms a single mag- 
netic domain (due to its shape anisotropy) and a sizable 
AB = (B L - B R ) /2 even if £ ext = [52]. This pre- 
viously unexplored regime proves to be highly beneficial 
for controlling DNSP. 

Even in the absence of on-chip magnets, the HFI be- 
tween thermally fluctuating nuclear spins and electrons 
creates an effective (Overhauser) field, B nuc , which sta- 
tistically varies between the two dots resulting in a small 
leakage current near B ext ~ and A ^ [33 , 34 . Com- 
pared to B nuc , the nanomagnet 's inhomogeneous field, 
B nm , is more stable in time, and |AB nm | ^> |AB nuc | (see 
Figure [1J). B ext is aligned along the easy axis (z-axis) of 
the nanomagnet, which has a coercive field of 52 mT. Be- 
cause the nanomagnet forms a single domain, B ext does 
not affect the magnitude of B nm . The associated AB nm 
results in a leakage current over a broad range of A and 
£? ext including distinct current maxima along the S di -T± 
resonances (Figure^). These current features are seen 
at small B ext and, hence, are not accessible with multido- 
main magnets (see above). Our observed current features 
are very different from measurements performed without 
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FIG. 2. (a) The relevant eigenenergies as a function of 
A, for Bnuc = 0. The inhomogeneous B nm causes singlet- 
triplet mixing (see enlarged S&-T- avoided crossing) and lifts 
the PSB. (b) Leakage current, /, through the DQD ver- 
sus A (stepped from top to bottom) and B ex t (swept at 
— 50mT/min). The red dashed line is a numerical prediction 
of the S a -T± resonances without DNSP using t c = 12ueV. 
The white dashed line includes DNSP (via the S a -T- res- 
onance) using eq [l] with r re i = 0.043 s _1 , 7 = 10 mT, and 
a = 35nA _1 s _1 . 



an on-chip magnet [33] . 

However, consideration of B nm alone does not explain 
all features in Figure [2}d. We must include the HFI and 
its effect on DNSP, as has proven necessary in other ex- 
periments pjJl I35H37] . The dashed red line in Figure [2]o 
is a prediction of the position of the S a -T± resonances as 
a function of £> ext and A. It takes into account B nm , but 
neglects B nuc [33 . Compared to this prediction, however, 
the measured resonances (blue arrows in Figure]^)) occur 
at larger |£> ext |. As we will show, this shift can be ex- 
plained by including DNSP, which produces a B* uc that 
compensates £? ext , e.g., B* uc < when £> ext > 0. 

The connection to DNSP becomes evident with the 
data shown in Figure [3^i, which probes the S a -T_ reso- 
nance as a function of time for a fixed A. We prepared the 
nuclear spin polarization, P, to P ~ by waiting three 
minutes at / = before turning on the voltage across the 
DQD. The current maximum, J max , at the S a -T_ reso- 
nance occurs later in time at larger B ext . Again, this can 
be explained if B* uc < and compensates B ext . Because 
the GaAs g-factor is negative, P = -5^/5™^, where 
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FIG. 3. (a) Current versus time, /(£), as a function of Bext 
along the horizontal line in Figure]^ (we prepared P(0) = 0). 
Density Plot: I(t,B ext ). The white dotted line predicts the 
moment of maximal current, t (/max), using eq[l]and the same 
parameters as in Figure J2jx Red Trace: I(t) measured at 
Bext = 25 mT. (b) P(P) for three different B ext using eq[l] 
and the same parameters as in Figure [2j). Far from resonance 
an exponential decay P — — r re iP remains (gray dashed line). 
#,•,• mark the B ext -dependent "adjustable" fixed point and 
■ a "trivial" fixed point, which appears near P = when 
#ext>45mT. 
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6.1 T is the Overhauser field magnitude pro- 
duced when all nuclear spins are aligned [33]. B* uc < 
implies that P > 0, which can only be explained if DNSP 
from T + outweighs that from T_ , despite the system be- 
ing near the S a -T- resonance. This peculiar situation 
results from spin-selective lifting of the PSB (of T_) and 
bolsters DNSP, as discussed below. 

Our explanation starts with a rate equation model p~9j 
[20l [38] including only polarization generated by T + near 
the S&-T- resonance (A comprehensive calculation in [33] 
Sec. Ill includes all (11) states). As a simplification, we 
use the average polarization P = (P L + P R ) /2. (P L + 
P R would mainly affect the decay of To, not T±.) The 
overall rate equation is 



Pit) = r pol (i) [i - p(t)} - r iel p(t), 



(i) 



where the polarization decay rate, T re i, is constant, while 
the build-up rate, T po i(t), is proportional to current, 
r po i(£) = &I(t) > 0, as observed experimentally. For 
convenience, we describe the current maximum at the 
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where J max is the (measured) resonant current and 7 is 
the effective width of the resonance. (Nonresonant states 
contribute weakly to I(t).) Here 



(3) 



is the energy of S a . We approximate E z by only including 
the average z-component of B 



inn? ^nm) so that 



E~ « \g\ MB (i? ext (t) - B™P(t) + B 



(4) 



Example P(P) curves are plotted in Figure [3)3 and are 
used to model the data in Figure [3^i. P (t = 0) = in this 
measurement, and the model predicts P(P = 0) > (evi- 
dent in Figure [3}}). Therefore, P increases in time until it 
reaches a stable fixed point at P = (and dP/dP < 0). 
For P ext < 43 mT, P passes through the S&-T- reso- 
nance, which coincides with P ma x 

in Figure u3b. As P ext 
is increased, the 5 a -T_ resonance moves to larger P, and 
with it move P max and the stable "adjustable" fixed point 
(A-FP, circles in Figure J3Jd). Accordingly, the measured 
(resonant) J max in Figures [3^l appears later in time with 
increasing B ext . When P ext « 45 mT, a second stable 
"trivial" fixed point (T-FP, square in Figure ^$jp) appears 
near P = and remains there for P ext > 45 mT. Hence, 
we expect P to remain near zero (far from resonance) 
at the T-FP. Indeed, no resonant current maximum is 
observed for P ext > 43 mT in Figure [3^i. 

Eq [I] provides quantitative predictions of the time evo- 
lution of the S&-T- resonance associated with the mea- 
sured /max- Namely, it yields the white fits in Figures^ 
and [3^i. These two separate fits share altogether four fit 
parameters. The agreement between our model and data 
indicates that the model captures the DNSP in both ex- 
periments. In addition, r re i agrees with reported values 

Our model reveals a straightforward procedure to max- 
imize P. We start at small P e xt where the T-FP is ab- 
sent and P is initialized at the A-FP (see top panel of 
Figure |4|l) . This initialization requires small P e xt and 
sufficient singlet-triplet mixing and is only possible with 
a single-domain nanomagnet. To reach a large P, we 
increase P e xt (with a sufficiently slow sweep rate) drag- 
ging the A-FP, and P along with it (see middle panel 
of Figure [4^). This dragging procedure works up to a 
maximum polarization, P max , occurring when the decay 
of P overwhelms its build-up and P m ax = (see bottom 
panel of Figure [4^) . P ma x is defined by solving eq [I] for 

P(Pmax)=0: 
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FIG. 4. (a) P(P) for three P ex t using eq[T]with 7 max = 1.2 pA 
and the fit parameters from Figure [2j). ■ marks the T-FP 
and # the A-FP. (b) Leakage current, /, for various A in the 
PSB regime [see (c)] measured while continuously increas- 
ing Bext (at 30mT/min) after initialization at the A-FP at 
Bext = 50 mT. (All but the lowest trace are vertically off- 
set in increments of 125 £A.) Sharp current peaks (marked by 
arrows) correspond to the maximum polarization, P ma x(A), 
in (a), (c) Pmax/ (1 — Pmax) from (b) (filled circles) and / 
(black/green traces), measured along the (black/green) lines 
in the inset, versus A. Inset: I measured as a function of 
gate voltages Vl,r (see Figure [Ir). The region of suppressed 
/ within the double triangle of finite current marks PSB. (d) 
Pmax/ (1 — Pmax) versus /, extracted from the traces in (c). 
The gray square is predicted from the data in Figure [2}). The 
black line expresses eqJHJwith a/r re i = 0.81 pA -1 using the 
fit parameters from Figure [2b. 
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Pmax > and both depend on A via I max . 
(the field corresponding to P ma x), the A-FP 
coincides with the S a -T_ resonance, and we expect a 
current maximum. 

Indeed, the leakage current measured during such 
DNSP sweeps, displayed for various A in Figure^, fea- 
tures a sharp current maximum at a particular P ext (see 
arrows), which we identify as P^ X (A). For P ext > 
P^ x , the A-FP is lost, and the nuclear spin polariza- 
tion decays (P < 0) with rate r re i. 

-Pmax (A) compensates P™^ X (A) so that we can equate 
P max (A) = B-r(A)/£C- • W(l - ^max) and / 
measured versus A are shown side-by-side in Figure ^ 
demonstrating a striking correlation between P max /(1 — 
Pmax) and I. In fact, in accordance with eq[5j Figure 
shows that P ma x/ (1 — Pmax) ^ I confirming our assump- 
tion that r po i ex I. The straight black line in Figure [4|i 
expresses eq [5] using the fit parameters from Figure |2p. 



The ability to explain three very different data sets (Fig- 
ures [2)3, [3^i, [4]i) with one set of fit parameters corrobo- 
rates the interpretation of the current peaks in Figure 
and the validity of our rate equation model [33] . 

Our highest P™? x — 2 - 9T (green data in Figure ^p) 
corresponds to P ~ 50 % (and generates an Overhauser 
field gradient of ~ lT/lOOnm across the DQD bound- 
ary). This exceeds by far previously reported polariza- 
tions in laterally defined DQDs [TQUT3] (a complementary 
measurement of P in [33 Sec. VI). 

To detail how B nm and the HFI combine to lift the 
PSB and induce DNSP, we compare our system with 
two simpler scenarios. If the HFI were the only mech- 
anism to lift PSB, no DNSP would be expected since 
all triplets are loaded equally often resulting in as many 
up as down nuclear spin flips. In experiments without a 
nanomagnet [36 , 39 , cotunneling weakly lifts the PSB (in 
competition with the HFI) and does so equally for each 
triplet, nearly irrespective of its energy. In contrast, the 
hyperfine-induced decay rate is strongly energy depen- 
dent. Therefore, near the 5 a -T_ resonance, T_ generates 
more nuclear spin flips than T + , and P < is observed 
without nanomagnet. 

In our case, B nm mixes T_ and S a strongly near their 
resonance, resulting in two (11) states that are no longer 
in PSB. Hyperfine-induced decay is heavily suppressed in 
these mixed states. In this situation, the HFI still con- 
tributes to the decay of T + (and T ) thereby producing 
DNSP and P > 0. In an alternative approach, DNSP 
has been studied for large t c (~ 100 ueV) by sweeping 
A [36 , [37j |39] . However, when t c ~ 1-10 peV, which is fa- 
vorable for spin qubits, P is limited by the energy of S a 
in the PSB regime, so that |P| < t c J (2 \g\ /i B £™ u a c x ) ~ 
10%. Moreover, B* uc < in our system provides a 
distinct advantage because B* uc compensates P ext such 
that the total effective magnetic field is constant during 
the polarization build-up; therefore, P is only limited by 
a^max/r re i when dragging P with B ext . 

We have demonstrated a nuclear spin polarization of 
^ 50 % in a DQD based on the enhanced ability to ma- 
nipulate the nuclear spin ensemble using an on-chip nano- 
magnet. Larger polarizations can be expected upon fur- 
ther optimization of the electronic spectrum, sample ge- 
ometry and materials. Our results demonstrate the flex- 
ibility offered by an on-chip nanomagnet, which could be 
used for all-electric ESR [29 while simultaneously polar- 
izing the nuclear spin ensemble at small t c values ideal 
for spin qubit operation. Such a system could be used 
to improve nuclear state preparation techniques [I2j l40l 
142] or for measuring complex nuclear phenomena such 
as spin squeezing [43 , quantum memory |23j [24] , dark 
states [44] , quantum phase transitions [45 , and superra- 
diance [46] . 
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I. OVERVIEW 

The following supplementary material provides additional information related to various aspects of the main article. 
We start in Section I^with details about the sample design and the experimental setup of the measurements. In Section 
11^] we introduce a model Hamilton operator which describes the hyperflne interaction (HFI) in our double quantum 
dot (DQD) setup including the inhomogeneous magnetic field of the nanomagnet. Based on a rate equation model, 
in Section IV[[ we perturbatively solve the dynamic nuclear spin polarization (DNSP) problem by explicitly taking 
into account the contributions of all four (11) states. We show that the perturbative solution justifies the simplified 
model used in the main article. Section V] provides detailed explanations of the current features in Fig. 2b of the 
main article. In Section V^]we present results of a complementary experiment that gives additional evidence of the 
validity of our model and of our interpretation of the data in Fig. 4 of the main article in terms of a large nuclear 
spin polarization. Section VI!] describes the fitting procedure for the data in Figs. 2b, 3a of the main article. 



II. SAMPLE DESIGN AND EXPERIMENTAL SETUP 



The samples have been fabricated from a GaAs / AlGaAs heterostructure containing a two-dimensional electron 
system (2DES) 85 nm below the surface. At cryogenic temperatures, the 2DES has a carrier density of 1.19 x 10 11 cm -2 
and a mobility of 0.36 x 10 6 cm 2 V -1 s _1 . Metallic gate electrodes (30 nm gold on top of 5nm titanium) have been 
fabricated on the sample surface by electron-beam lithography and standard evaporation/lift-off techniques (Fig. SI i) . 
The Co nanomagnet with a thickness of 50 nm was evaporated directly on top of the leftmost gate and capped with 
5 nm of Au to prevent oxidization. Negative voltages applied to these electrodes are used to deplete locally the 2DES 
and thereby define the DQD. The absolute electron occupation, (m n), was determined by quantum-point-contact 
charge detection |Rlj . All measurements have been performed in a dilution refrigerator at an electron temperature of 
~ 100 mK. Fig. SI 3 sketches the experimental situation in this Letter. A source-drain voltage of V = (/is — Md)/c 
is applied across the DQD between degenerate leads. The DQD is in the Pauli-spin blockade, where a triplet state 
can only contribute to current if it is coupled to a singlet state, e.g., by field inhomogeneity or interaction with the 
ensemble of A/" L ( R ) ~ 10 6 nuclei. 



III. THE HAMILTONIAN 



The total Hamiltonian of the system includes electrostatic, magnetic, and hyperflne contributions and is given (in 
the relevant subspace depicted in Fig. [sT] o) by 

H = H e \ + Hb + H hi . (SI) 
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FIG. SI. Experimental setup a) SEM image of the DQD device used in the experiments, b) A DQD in Pauli-spin blockade 
(typical experimental situation) , while a voltage V is applied between the degenerate 2D leads with chemical potentials /is and 
/xd- Vertical lines are tunnel barriers. The right dot always contains at least one electron. The dashed horizontal lines are the 
spin-split chemical potentials of the spin up/down (0 1) states, where the chemical potential of a quantum dot is defined as the 
energy needed to add one more electron. The solid horizontal lines are the chemical potentials of the five relevant two-electron 
basis states, where Ez is the Zeeman energy and A is the energy detuning between the singlet states, Sn and So2- 



Using the diabatic singlet-triplet basis {T + , T , T_, 5n, 602}, the electrostatic contribution is 

H el = t c /2 (\S n ) (S02I + |5 02 > (S11I) - A |5 02 > (S02I , (S2) 

where A is the interdot energy detuning (see Fig. |Sl) ) and t c is the interdot tunnel splitting (see Fig. 2a of the main 
article). The interaction between the local magnetic fields and the electron spins in the two dots is described by 



H B = gfJLB [(Bext + B^ m ) • S L + (B e 



(S3) 



where B ext is the external magnetic field, the local magnetic field of the nanomagnet, S L ' R the local electron 

spin operator, g the electron g-factor, and [1b the Bohr magneton. 

The wave function of an electron confined in a lateral GaAs-dot overlaps with ~ 10 6 nuclei. The hyperfine interaction 
between the electron spin and the nuclear spins is dominated by the contact term i^ R = A^ ,R I^ ,R • S L,R , where 
I^' R is the fcth nuclear spin operator and A^' R the hyperfine coupling constants. A^ ,R is proportional to the overlap 
between the wavefunctions of the kth nucleus in left/right dot and the electron and varies by isotope type. It is 
common to define an average A reflecting the natural abundance of isotope type and the average overlap with the 
electron wavefunction. In this approximation, the contact Hamiltonian is i^ R = AI L ' R • S L,R , where I L ' R is the 
average nuclear spin (ensemble) operator and A = 85ueV in GaAs[R2 . Electrons in the left /right dot couple to 
different sets of nuclei, and we can write 



Hm — H 



Hi 



A (I L • S L + I R 



i=L,R 



PSl 



p + si + p_s\ 



S^ R ± iS^ R and 7±' R 



(S4) 



ib R ±ub R 



where 5^' R and Z^' R are the corresponding z-projection operators; and 5^' R 

are the spin raising and lowering operators. In a semiclassical approximation I L,R can be replaced by the effective 
nuclear magnetic (Overhauser) field |R3j B^ R = A (l L ' R ) / (^b), where (...) denotes the expectation value, and 
(7) max = 3/2 in GaAs. In ESR experiments (not shown), we measured g ~ —0.36 in our system. This predicts, 
for fully polarized nuclear spins (P = 1), an Overhauser field magnitude of B™* = A (/) max / (\g\ /jLb) — 6.1 T. The 



semiclassical version of has the same form as Hb (see Eq. (S3)), and we can summarize 

H B + ffhf = 5MB (B L • S L 



(S5) 



where B L ' R = B ext + B^ + B^ R is the total effective magnetic field acting on an electron in the left and right dot, 
respectively. 

In analogy to B = (B L + B R ) /2 and AB = (B L — B R ) /2, we define the symmetric and antisymmetric spin 
operators S = (S L + S R ) /2 and AS = (S L - S R ) /2. We then use B± = B x ±iB y , AB± = AB x ±iAB y , S± = S x ±iS y 
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and AS± = AS X ± iAS v , denned akin to the spin raising and lowering operators in equation (S4), to write equation 



(S5) analogous to the right hand side of equation (S4): 

H B + H M = gi i B (2B Z S Z + 2AB Z AS Z + AB + AS_ + AB_ AS+ + B+S- + B-S+) . 



(S6) 



With the quantization axis, z, defined parallel to B ext , the matrix representation of the (semiclassical) total Hamil- 
tonian in the basis spanned by the diabatic singlet and triplet states {T + , T , T_, 5n, S02} is 
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(S7) 



where n* = t* = t c /fjL* and A* = A//i*. The matrix representation (S7) illustrates that the x- and y- 

components of AB mix T± with 5n, while the z-component of AB mixes Tq (which has no spin component along 
the z-axis) with Su. Instead, B z leads to the Zeeman splitting of the T± states. Note that the off-diagonal terms 
B±, which mix T± with Tq, vanish if the quantization axis is chosen parallel to B. 



IV. PERTURBATION CALCULATION OF DNSP 



In the main article, we use a simplified approximation that only considers the hyperfine contribution from T + . In 
this section, we calculate the hyperfine- induced decay of all (11) states using perturbation theory (Fermi's golden 
rule) similar to reference |R4j in which, however, the effects of a nanomagnet were not included. Here we show that 
the simplified model produces the pertinent features of the perturbation theory, justifying the approximation used in 
the main article. 



We start by writing equation (SI ) as H = Hq + + H^ , where 

H = H el + H B + 2gp B (XucSz + AB£ UC AS*) , (S8) 
Hi = gf i B (Bt c S T + AB± UC AS T ) , (S9) 

are the bare Hamiltonian and hyperfine flip-flop Hamiltonians, respectively. We treat H^ + as a perturbation 
of Hq. Diagonalization of Hq provides the unperturbed eigenvalues, E n , of the nth eigenstate, \n). We account for 
coupling to the leads by a simple master equation with four Lindblad operators (eliminating the intermediate stage 
in the sequential tunneling process (0 2) — >• (0 1) — >• (1 1)) and assuming that the four (1 1) states {T + , T , T_, Su} are 
populated with equal rate: 

j t p = ±_[H 0: p] + ^ (\ x ) ( S 02\p\So2) (x\ - \ (|S 2> (S02\P + P\S 2) <5 02 |)) . (S10) 

We approximate the dynamics by a rate equation for the populations p = (pn, P55) m the five energy 

eigenstates. 

f = %. (Sn, 

The transition matrix = (G^)ij, describes decay of the level n with a rate = r#s n , determined by s n = 
I (£02 1 n) | 2 , the overlap of \n) with the localized singlet state. Since only (1 1) states are refilled, the rate with which 
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\m) is populated is proportional to (1 — s m ): Hence the matrix elements of are 

G ( °l= ^f(l-s n )s m (n^m), (S12) 
Gi°2 = - I fs n (3 + Sn ). (S13) 

The width of the level \n) is given by HT^ = HTjis n . 

We include the hyperfme flip-flop processes using Fermi's golden rule (assuming a constant density of states over 
the range of energies of the (11) states) to determine the flip-flop rate from an initial state, |n), to the final state, |/) 
as 

where the factor (1 =p P) /2 expresses the influence of polarization on the nuclear-spin flip rates. The total escape rate 
from \n) is then 

rn = r° + X)(r+_ / + r-_ / ). (su) 

We can neglect cotunneling in our sample, since its contribution to the leakage current is negligible compared to that 
of the nanomagnet. The full rate equation is now given by 

f= G p, (S15) 

with G nm = G ( nl + T+^ m + T-^ m (n ^ m) and G nn = g1°2 - E m / n ( r n-^m + r n^m)- 

We determine the leakage current and the DNSP rates by numerically solving for p = 0, obtaining the steady state 
populations p^. The magnitude of the total current is given by 

I = I°+I+ + I- = e J2pnVn, 

n 

where e is the magnitude of the electron charge. 1° is the non-polarizing current induced by the nanomagnet, while 
J ± are the hyperfine-generated currents that polarize in opposite directions. 7 ± are expressed as 

I ± =eJ2PuTt. 

n 

To convert these currents into nuclear polarization rates, we write 

Pol eN ' 

where we have normalized by e and included that the polarization, P = (7V^ — N^) /TV, changes by 2/N per nuclear 
spin flip for TV nuclei in total. The overall polarization rate equation is 

> = r+ 1 -r- 1 -pr rel . (sie) 

In the main article, we neglect r~ ol and approximate 

P - Tpol - ^rel - Otl (1 - P) ~ PT Ieh (SU) 

where a is taken as a constant (~ 35nA _1 s _1 ). 
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To compare the above theoretical model with the simplified version used in the main text, Fig. S2 shows P versus 
P and B^t predicted using the two models represented by equations (S16) and (S17). The simple model, shown 



in Fig. |S2p ,, predicts a narrow region of positive P running diagonally through the P-5 ext plane. The perturbation 
theory calculation, shown in Fig. |S2| p, predicts this same feature, which is associated with the S a -T_ resonance, as 
well as more complicated behavior associated with other singlet-triplet resonances (see Fig.|S2fc). Both models predict 
that the S a -T- resonance is the only resonance which creates a stable fixed point at large P and B ext . Figs. |S2]>e 
demonstrate the good agreement between the simplified model (orange lines) and the perturbation calculation (black 
lines) at the /S a -T_ resonance. (The orange lines here are the same as the lines in Fig. 4a of the main article.) 
Compared to the other two resonances, P at the S a -T± resonances is very steep and provides strong feedback of the 



nuclear spins toward the A-FP, a distinct advantage for DNSP. From Fig. S2, it is evident that equations (1-4) of the 
main article are sufficient to model the nuclear polarization associated with the S a -T- resonance. 

The leftmost and rightmost P extrema in Figs. S2 >d are generated by the S s -T± resonances, while the small feature 
between the S a -T± resonances corresponds to the crossing of the triplet levels at B z = B ext + B nm — PB™* « 0. 
Note that the nuclear field components B^ C ,AB^ C ,Z play an important role for our treatment, determining directly 
the strength of the hyperfine flip-flop rates. For the DNSP rates depicted in Fig. [S2] we have averaged calculations 
for different values of the nuclear field fluctuations chosen according to a Gaussian distribution with zero mean and 
standard deviation of ~ 3 mT. 
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FIG. S2. Nuclear polarization rate a) The nuclear polarization rate, P, calculated analytically using equation (1) of the 
main article. The dashed horizontal lines correspond to the constant B eyLt slices shown in Fig. 3b of the main article, b) P 
calculated numerically using equation (| S16| with Tr = 4.5 GHz and N — 10 6 nuclei. A random Gaussian distribution was 
used for B^ uc and B^ uc as in reference |R5| . The scale bar applies to both a and b. Both calculations were performed using 
t c — 12ueV, A = 8yieV, r re i = 0.04 s _1 , and the magnetic field distribution of the nanomagnet (Fig. lc, main article), c-d) P 
as a function of P at constant B ex t (lines in a,b) co mpa ring the rate equation model [orange, equation (1) of the main article] 
with the perturbation calculation [black, equation (S16)]. The ranges in P of these slices are indicated by the magenta lines 
in b. The polarization rate extrema associated with the four singlet-triplet resonances are labeled in c. The top axis indicates 
the z-component of the total effective magnetic field, B z = B e ^ t + B nm — PB™*. In addition to the extrema associated with 
the singlet-triplet resonances small polarization features appear near B z — (dashed vertical lines) where triplets become 
degenerate. 
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V. HYPERFINE- AND NON-HYPERFINE-INDUCED LEAKAGE CURRENT 



The leakage current through the DQD shown in Fig. 2b of the main article contains a number of features which 
can be traced back to the inhomogeneous field produced by the single domain nanomagnet. To illustrate this, in Fig. 
S3 we show the PSB leakage current through two different DQD devices with identical gate layout. The data in Figs. 
S3 



i and [S3p show measurements for opposite sweep directions acquired on the sample also presented in the main 
article but for a larger interdot tunnel coupling t c ~ 15ueV. These data are more richly featured compared to those 
in Fig. S3 z in which no nanomagnet was present (t c ~ 1 ueV). 



In Fig. |S3p , the measured current is approximately symmetric with respect to the B ext = axis, and the main 
features are increased current along the A = and B ext = axes (for A > 0, that is, outside of Coulomb blockade) 
and a global maximum at A = £? ext c± 0. Similar data have already been published and discussed in detail in reference 
|R6j . In short, current is created by t c in combination with the hyperfine interaction, which mixes triplet and singlet 
states strongly when A ~ or £? ext — 0. The width of the current maximum along the B ext = axis is determined 



by the standard deviation of the fluctuating B nuc |R5l IR6j . The data in Figs. S3i and S3d illustrate that the sizeable 
B nm adds complexity. The following list provides a short explanation for each of the features specific for the sample 
with nanomagnet: 



The most obvious response to expect when sweeping B ext is hysteresis of the magnetization of the single domain 
nanomagnet. Because of its single-domain character, we expect the nanomagnet to switch polarization abruptly 
when B ext passes the coercive field. An abrupt switch in magnetic field, in turn, should cause an equally abrupt 



change in the leakage current signal. Such features are indeed observed in Figs. [S3p,b at £? ext ~ ^52 mT, 
respectively (see white arrows), and are also seen in Fig. 2b of the main article. 

In the presence of B nm , the eigenenergies of the T±-like states are never zero. However, the relevant magnetic 
field, |B L + B R | can be minimized by B ext , and at this minimum, the T± are most degenerate and a local 
maximum of the leakage current is expected. For our B nm values, |B L + B R | is minimized at B ext ~ ±12 mT, 
depending on the polarization of the nanomagnet. These fields are indicated in Figs. |S3^i,b with red arrows and 
faithfully identify the current maxima. 



The observation of distinct local current maxima at the S Ei -T± resonance (black arrows in Figs. S3i and[S3p) 
and in Fig. 2b of the main article) is unique to samples containing a single domain nanomagnet. The sharpness 
of these peaks can only be explained by taking into account the hyperfine induced dynamics of the nuclear spins. 
The actual position of the S a -T± resonance is shifted towards larger |i? ex t| compared to its prediction (see Fig. 
2b of the main article). In the main article we explain this shift by taking into account the hyperfine induced 
DNSP. 



mixes triplet and singlet states weakening the spin blockade and allowing leakage current to flow. However, 
B ext tunes this mixing. In fact, the condition AB || B defines a local minimum of the singlet mixture with 



the T± states |R5j . This can be readily seen from the Hamiltonian in equation (S7) if the quantization axis 



is chosen parallel to B. For our system, this condition is satisfied when £? ext ~ ±8mT. We actually observe 



current minima at slightly shifted values (see yellow arrows in Fig. S3i,b) owing to the complex DNSP that 



occurs while sweeping £? ext (see Fig. S2). 



The current at finite B ext and small A in Figs.[S3p and[S3p is characterized by strong switching noise and dragging 
effects which has also been observed in samples without on-chip magnet [R61IR7] . We for so a detailed discussion of 
these effects, which can be explained in terms of the hyperfine dynamics in the presence of more than one stable fixed 
point at A - 0[R4], 
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FIG. S3. Spin-blockaded leakage current, a) The dc leakage current, /, as a function of B ext (swept from positive to 
negative) and A (stepped from positive to negative) (t c ~ 15 ueV. b) Same as in a but for the opposite sweep direction of B ex t 
(from negative to positive). Arrows in a and b mark specific features explained in the main text, c) / as a function of B e *t 
and A measured using a sample with an identical gate layout as the sample used in the main article, but without the on-chip 
nanomagnet. Overall, maximum / values are one order of magnitude larger than in a and b owing to stronger source/drain 
coupling with this particular gate tunning, but the region of enhanced / is much smaller. The magnetic field was stepped from 
right to left and the energy detuning was swept from bottom to top. The perpendicular sweep direction compared to a and 
b does not affect the main features of this measurement. It does however cause a noisy background which is typical for this 
sweep direction and is caused by charge noise triggered by sweeping gate voltages. 
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VI. NUCLEAR POLARIZATION BUILD-UP 



Here we present additional data to support our interpretation of the DNSP data. Fig. [S4] demonstrates that the 



polarization dragging in Fig. 4b of the main article is reproducible. All the main features in Fig.[S4fi are reproducible, 
especially the position of B™* where polarization is lost. The current measured at the beginning of each field sweep 
near £? ext = is the typical leakage current that appears near £? ext = (see, for example, Fig. S5) and is extended 



somewhat because of DNSP. We interpret the sharp current maximum (labeled B™*) as the point of maximum 
polarization (see main article) . The SI in Fig. S4 3 results from losing the stable polarization condition and related 



resonant current as the polarization decays and the system drifts away from resonance. 



Fig. [55] shows a second technique for demonstrating the polarization created during fixed-point dragging measure- 
ments. In Fig. [S5^ l, we show I measured as a function of B ext and A over a much larger range of B ext than Fig. 2b 
of the main article. The current features at low B ext in S5i differ from those in Fig. 2b of the main article mostly 
because A was swept rather than B ext . 

The current traces in Fig. [S5} d, measured at B ext = 250 mT after the nuclei have been polarized to P ^ 

0%. 



= 17.5 mT and P 



4%, 
The 



are very similar to the trace in Fig. S5 3, which shows current measured at B ext 
quantitative similarity between the P ~ 4% current traces in Fig. |S5]d and the trace in Fig. S5^ allows us to conclude 
that B* uc compensates B ext (B* uc ~ —B ext for B ext ^> |B nm |) thereby reducing the total effective field. In contrast 
to the polarized traces in Fig. [S5p, which are almost symmetric with respect to A = 0, the P = % curve in Fig. 



S5z is asymmetric and exhibits switching noise for A > 0. We attribute this behavior to small changes in nuclear 



spin polarization, while in Fig. S5 3, the polarization is stabilized at the adjustable fixed point (A-FP). The current 
traces in Fig. |S5| o measured at P « 4% are repeatable and much larger than the trace measured at P — 0%. This 
demonstrates that the polarization is finite and stable. 



Current features, such as the local maxima in Fig. 2a of the main article and Fig. S4 at B ext > 1 T, are not seen 



in Fig. S5i. These features are missing because A is swept, and significant polarizations are not obtained. Taken 



together with Figs. 2,3 of the main article, Figs. [S4| and |S5] demonstrate the ability of our system to generate large 
nuclear spin polarization — and detect it. 



500 



250 




FIG. S4. Polarization sweep repeatability a) Five polarization sweeps measured at a rate B e ^t — 35mT/min and 
A = 150 ueV. The traces are offset by N x 120 fA; N = 0, 1,2,3,4 for clarity. The current bistabilities observed beyond £™? x 
are consistent with DNSP[R4, R6 . b) The details of a current trace near B^ x taken from within the boxed region of a 
demonstrate a sharp resonance and clear change in current, SI c± 11 fA, before and after sweeping through 5™t x - 
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FIG. S5. PSB leakage current at large B e xt a) Leakage current, /, as a function of A and a large range of £> ex t- The 
magnetic field has been stepped and A has been swept from negative to positive to minimize DNSP effects. Inset: Enlarged 
region centered at Bext = 250 mT showing the small (~ 100 fA) leakage current when P = 0%. b) The solid orange lines are / 
measured while sweeping A multiple times through A = into Coulomb blockade along the orange arrow in a. These data were 
measured after ramping B e *t from Bext — to B ex t — 250 mT (along the white arrow in a) creating a polarization of P ~ 4%. 
The dashed blue line is / at P = % extracted from a at B ext = 250 mT and is negligible compared to / with polarization, c) 
/ measured versus A starting with unpolarized nuclei extracted from a at B e *t = 17.5 mT. 



VII. DATA FITTING 



The fitting procedure of the data in article Figs. 2b and 3a involves solving numerically the nonlinear differential 
equation given by equation (1) of the main article for a given set of parameters. This produces P(i), which is then 
fed into equation (2) to find I(t). Our goal is not to reproduce the details of the measured I(t) traces, but only the 
position of its maximum at the S a -T- resonance, that is, the position of the resonant current 7 max . Therefore, the 
final step is to calculate the position of / max using the numerical I(t). This procedure was repeated with different 
parameter sets until an agreement between theory and data was found. 

The numerical fit needs the following parameters: t c , 7, a, 7 max , and r re i (see equations (1-3) of the main article). 
I max and the ratio a/r re i = 0.8 pA -1 (see main article Fig. 4d) were measured, thus reducing the overall number of 
fit parameters to three. 

For the time-dependent data (see main article Figs. 3), 7 max = 100 fA is the measured peak height. i ma x values for 
the A-dependent data (see main article Fig. 2b) are unique for each value of A because I is A dependent. Fig. [S6] 
details how I max (A) is extracted from the / measured as a function of B ext and A. The main result is that I max (A) 
is identical to 1(A) measured near B ext = 0, where spin blockade is lifted. 

When B nm and B nuc are known, the A-£? ext position of the S a -T- resonance can be approximated analyt- 
ically. This approximation is used in equation (4) of the main article and includes only B* m providing ~ 
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we compare exact numerical results with the analytical approximations 
for all four singlet-triplet resonances. The analytical approximation for the S a -T- resonance is in excellent agreement 
with the numerical calculation for £? ext > 0. 

One set of fit parameters (t c = 12ueV, r re i = 0.043s _1 , 7 = 10 ± 1 mT, and a = 35nA _1 s _1 ) reproduces the 
data of two very different experiments in Figs. 2b, 3a of the main text. The data sets were measured using identical 
gate voltages. For a slightly different system tuning, only two parameters are expected to change, namely t c and 
a, because they reflect the various hyperfine and non-hyperfine system rates, which are strongly gate dependent. 7 
depends mostly on B nm , which is constant, while r re i should be independent of gate tuning because it is a property of 
the nuclei. These expectations are supported in Fig. [S8] where I versus B ext and time has been measured after making 
the voltage of the top center gate more positive (see the gate design in Fig. la of the main article). The ability to 
describe disparate sets of data in different tuning regimes with either no change or only justifiable adjustments to fit 
parameters demonstrates the validity of our model. 
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FIG. S6. Determination of 7 ma x. a) Leakage current, /, measured versus B ex t and A repeated from Fig. 2b of the main 
article. After searching within the boxed region, the • indicate the Bext-A position of / max along the S a -T- resonance, b) 
Here • are the 7 ma x values from a plotted as a function of A. In comparison, the black line is / measured versus A along 
the symmetry axis (black arrow in a). As expected, /max (A) along the S&-T- resonance follows the general /(A) trend. A 
smoothed version of /(A) is used to create the smooth numerical fit in Fig. 2b of the main article. 




FIG. S7. Location of singlet— triplet resonances. The numerically calculated data points indicate the position in the 
Bext-A plane where singlet and triplet states are resonant. The solid lines are approximated values calculated analytically. 
Each line is labeled to identify which states are resonant. In particular, the black line is the analytic approximation for the 
S&-T- resonance and was used to fit the DNSP data (see Fig. 2 of the main article). The exact and approximate values 
are in excellent agreement for B e ^t > 0, where all DNSP measurements were performed. The upper branch of the numerical 
calculation is included in Fig. 2b of the main article. The vertical dashed line at B e *t — — 12 mT, which is the minimizing value 
of |B L + B R |, defines the symmetry axis where the triplets are most degenerate. Here t c = 12ueV. 
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FIG. S8. Nuclear Spin Dynamics. Current measured as a function of time demonstrating DNSP at the S&-T- resonance. 
The dashed orange line is a numerical fit of the data. Here the system is tuned to a larger tunnel coupling compared to the 
tuning used in Fig. 3 of the main article. As a result, t c = 22ueV is larger, which makes E a larger, and the current peak is 
not observed until B e *t ~ 20 mT. Meanwhile, a = 14nA _1 s _1 suffers from the increased t c . The other fit parameters are 
unchanged, that is, 7 = 10 mT and r re i = 0.04 s -1 . 
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